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ABSTRACT 

In the present work we show a comprehensive comparison of our new stellar pop- 
ulation synthesis MIUSCAT models with photometric data of globular clusters and 
early-type galaxies. The models compare remarkably well with the colours of Milky 
Way globular clusters in the optical range. Likewise, the colours of M31 globular clus- 
ters can also be explained by the models by assuming younger ages then their Galactic 
counterparts. When compared with quiescent galaxies we reproduce the colour evo- 
lution at intermediate redshift. On the other hand we find that the colour relations 
of nearby early-type galaxies are still a challenge for present-day stellar population 
synthesis models. We investigate a number of possible explanations and establish the 
importance of a— enhanced models to bring down the discrepancy with observations. 

Key words: galaxies: elliptical and lenticular, cD — galaxies: stellar content - 
globular clusters: general — galaxies: star clusters: individual: M31 



1 INTRODUCTION 

Stellar population synthesis (SPS) models are com- 
monly used for solving a great variety of studies, from 
age/metallicity estimates in globular clusters to the recon- 
struction of the star formation history (SFH) in galaxies. 
The modeling of the spectral energy distribution (SED) 
emitted by evolving stellar populations requires at least 
three main ingredients: stellar evolutionary tracks, a spectral 
or photometric library with well-established atmospheric pa- 
rameters for each library star, and an initial mass function 
(IMF). A variety of models exists in the literature t hat differ 



in the i r treatment of e a ch of these ingredients, e.g . , IWorthev 
Vazdekis! 1 19991: jBruzual fc Chariot! 120051 : iMaraston 



199-1 



2005 



Conrov fc Gunnl l201fj ). In particular, the adoption 



of the spectral library used to convert physical quantities 
in observed flux turned out to be critical in the compari- 
son with observations. Although theoretical libraries (e.g., 
iLeieune. Cuisinier. fc Bused 1 19981) are widely used anion; 



the most popu lar SPS models (jBruzual fc Chariot! 120031 : 
iMarastonl l2005h because they are easier to handle, given 
the wide stellar parameter and spectral coverages and the 
arbitrary wavelength resolution, empirical-based libraries 
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represent an advantage as the stars are real. The adop- 
tion of empirical spectra allows to overcome the uncer- 
tainties in the underlying model atmosphere calculations, 
giving more r eliable predictions when com pared with both 
spectroscopic (jVazdekis et al.ll2003l. |2010|'1 and photometric 
IMaraston et al.ll2009l : |Peacock et al]|201ll ) data. 

In Vazdekis et al (2012, hereafter Paper I) a li- 
brary of single-burst stellar population model SEDs cov- 
ering the spectral range 3500 - 9500 A at moderately 
high resolution is presented. For building-up these mod- 
els we combine three exte nsive empirical stellar spectra l 
libraries, namely MILES ijSanchez-Blazquez et al] [2006). 
Indo-US dValdes et alJ2004l and CaT JCenarro et alj2001al . 
ICenarro et al.l l2001bl . Icenarro et al]|2007l ). Here we focus 
on calibrating these models by comparing these new predic- 
tions to a variety of photometric data ranging from globu- 
lar clusters (GC) to quiescent galaxies. Broadband colours 
are among the simplest predictions of population synthesis 
models, thus comparing the model predictions to observed 
colours is the natural zeroth-order test of compatibility be- 
tween models and observations. 

Since GC are almost coeval and mono-metallic popula- 
tion of stars that closely resemble a single stellar population 
(SSP), they are the ideal benchmark to test the accuracy 
of the models. On the other hand, quiescent galaxies can 
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reach larger metallicities than the metallicity distribution 
of most globular cluster systems and allow us, with some 
caveats, to calibrate the models even in the high metallicity 
regime. Although they can be modelled as an SSP as a first 
approximation, pieces of evidence of recent star formation 
in ea rly-type galaxies have been found by different authors 
(e.g.. [Bressan et alj|200d : ISarzi et al.ll2008h . which indicate 
the need of more complex stellar populations to obtain bet- 
ter fits. Moreover, galaxies may contain a certain amount of 
internal reddening that can further complicate the analysis. 

The paper is organized as follows. In Section 2 we 
present a brief summary of the stellar population synthe- 
sis models presented in Paper I. In Section 3 we describe 
the webtool developed to download and handle these mod- 
els. In Sections 4 and 5 we present the comparison of these 
models with globular cluster and galaxy data, respectively. 
Section 7 summarizes our results. Throughout the paper, we 
adopt the concordance cosmology: Ho— 70 km s^ 1 Mpc , 
m = 0.3 and Ov = 0.7. 



2 MIUSCAT MODELS 

The stellar population models we use throughout the pa- 
per have been presented in Paper I, here we briefly sum- 
marize them. These models represent an extension of the 
IVazdekis et al. I (|200l l2010h models, based on the CaT 
and MILES empirical stellar spectral libraries. MILES 
includes 985 stars covering the range AA 3525-7500 A 
( Sanchcz-Blazq uez et al.| [2006) with a spectral resolution of 
2.5 AFWHM, (Falcon- Barroso et al. 2011). The CaT library 
(|Cenarro et al.1 l2001al ) consists of 706 stars in the range 
AA 8350-90 20 A at resolution 1.5 A. The Indo-U.S. stel- 
lar library (|Valdes et al.|[20o3 ) has also been added to fill- 
in the spectral range around 8000 A that is not covered by 
the MILES and CaT libraries and to extend blueward and 
redward the wavelength coverage of the models. The Indo- 
U.S. includes « 1200 stars with spectra that typically cover 
from 3465 to 9469 A at resolution ~1.36 A FWHM (Falcon- 
Barroso et al. 2011). 

Although all the stars of the Indo-U.S. library have been 
used to obtain an homogenized set of stellar parameters 
matching those of MILES and CaT libraries, only a sub- 
sample composed of 432 stars have been included to build- 
up the final MIUSCAT catalogue to feed the models. Only 
stars without gaps in the stellar spectra or significant telluric 
absorption residuals, among other reasons, were considered 
(see paper I for details). As the flux-calibration quality of 
the Indo-U.S. is not as accurate as in the other two libraries, 
a procedure was applied to correct the spectrum shape in 
the region AA « 7400-8350 A to match the MILES and CaT 
stellar spectra. Finally the Indo-U.S. spectrum was used to 
extend the spectral range blueward the MILES and redward 
the CaT spectral ranges without applying any correction 
to the Indo-U.S. spectrum. Finally, the CaT and Indo-U.S. 
spectra were smoothed to match the resolution of MILES, 
i.e. 2.5 A (FWHM). 

The MIUSCAT st ellar spectra were implement ed in the 
models as described in IVazdekis et alj |2003l . l20ld ) and up- 
dated in Falcon-Barroso et al. (2011). Finally the relevant 
spectral ranges of the resulting SSP SEDs were extracted 
and plugged to the original model SEDs computed on the 



basis of the MILES and CaT libraries. For the latter step 
a further correction, similar to that described above for the 
individual stellar spectra, was applied to the model SEDs to 
obtain a unique spectrum, which is identical to the MILES 
and CaT based models in the corresponding spectral ranges, 
respectively. 



2.1 U magnitude 

Unlike MILES SSP spectra, the new MIUSCAT models al- 
low us to measure from B to I broad-band filters and related 
colours. It is als o possible to derive th e U filter, which starts 
at » 3050 A l|Buser fc Kuruczlll978h . by correcting for the 
missing flux blueward 3464.9 A. For this purpose we make 
use of the fractions provided in paper I, which are obtained 
with the aid of the library of stellar spectra of |Pickleslll998l 
(AA 1150-25000 A). For each SSP the fractions obtained for 
the different stellar spectral types are integrated along the 
isochrone to obtain the corresponding missing flux for the 
SSP SED. Finally the total U magnitude is derived by mea- 
suring the flux of the MIUSCAT SSP spectrum through the 
U filter redward 3464.9 A and correcting by the missing 
flux. To compare our model predictions to the observational 
data discussed in this paper we follow the same method to 
compute the SDSS u filter at two redshift values (i.e. and 
0.04; see Section 6). 



3 THE WEBTOOL 

The model predictions used in this work have been obtained 
with the aid of the webtool facilities that we provide in the 
MILES website: |http://miles.iac.es| The MILES and CaT 
stellar libraries and the stellar population synthesis mod- 
els can be retrieved and handled according to the user re- 
quirements. The webtools include the transformation of the 
spectra to match the instrumental set-up of the observa- 
tions (spectral resolution and sampling), measurements of 
line-strength indices and synthetic magnitudes derived from 
the spectrum. 

As in this work we apply the MIUSCAT models to a 
large sample of galaxy data we describe in more detail an- 
other webtool facility that will be extensively used in Sec- 
tion 6. The tool returns the model spectrum corresponding 
to a given Star Formation History (SFH), which can be ei- 
ther user-defined or follow some parametric description. The 
first case is devoted to handle with complex SFHs, with an 
arbitrary number of bursts. For instance, this case will be 
particularly useful for the SFHs retrieved from cosmologi- 
cal simulations. Given the mass fraction and the metallicity 
at different ages, the tool computes the corresponding SED 
according to: 



m,iF\(U, Zi) 



(1) 



where Nb is the number of input bursts, m, is the fraction 
of mass formed in the i-th bin and F\(ti,Zi) is the i-th 
SSP spectrum for a chosen IMF. In the case that ti and Zi 
are not included in the age /metallicity grid of the models, 
the closest SSP model to the desired parameters is selected. 
Thus no interpolation is performed on the SSP spectra. 
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For the parametric SFHs four different cases are al- 
lowed: 

• Multiple bursts: up to five bursts occurring at different 
ages and with different metallicities are summed up 

• Truncated: the star formation rate (SFR) is assumed 
constant between the time of formation, tp, and the trunca- 
tion time, tr (all the times are referred to the present time) 

• Exponential: the SFH is modelled as a power-law rising 
term at early times plus an exponentially declining SFR: 

SFR(t) = (!£_£)»e— (2) 
t F 

where tF is the time when star formation starts, r is the 
e-folding time and n parametrizes the power-law rising term 
at early times. 

• Exponential plus bursts: up to five bursts can be added 
on top on the exponential SFR 

As we will show in the following sections, the SFH 
webtool is particularly useful to construct galaxy SEDs, for 
which the use of complex stellar populations provide better 
fits than the SSP approach. In Section 6 several applications 
of the parametric SFH will be explored. 



4 COMPARISON WITH GLOBULAR 
CLUSTERS 

4.1 Milky Way globular clusters 

Here we use the catalogue of iHarris! Il996l (2010 edition) 
as a reference for the Milky Way (MW) globular clusters. 
The catalogue is a compilation of 150 clusters including 
the Johnson-Cousins UBVR1 photometry, [Fe/H] metallicity 
and reddening E(B— V) among other data. Colours, except 
B— V, have been de-reddened usin g the extinction law of 
ICardelli. Clayton, fc Mathisl (|l989h . 

Figure [1] shows the comparison of the optical broad- 
band colours of the GCs with the ones derived from the 
models. The model metallicity indicates the total abundance 
ratio [Z/H]. A fit to the data points with a second order 
polynomial is indicated to guide the eye. In order to avoid 
clusters with uncertain photometry due to high extinction, 
we excluded the most reddened clusters with E(B— V) >0.5 
(open circles). The M1USCAT models are shown for ages of 
10.0 and 12.6 Gyr for a Kroupa Universal IMF. The use of 
a Salpeter IMF would leave the colours almost unaffected 
(with differences smaller th an 0.01 mag). To derive the syn- 
thetic colours we adopt the lBuser fc Kurucj |l97^) filter re- 
sponses and Vega magnitu des. The zero point has been set 
by using the iHavesl ll 19851 ) spectrum of Vega with a flux of 
3.44 x 10~ 9 erg cm" 2 s" 1 A -1 at 5556A(see Falcon- Barroso 
et al. 2011 for further details on the computation of synthetic 
magnitudes) . The colours of Vega in the Johnson-Cousins fil- 
ters are assumed to be U— B=U— V=B— V=V— R=0 mag, 
V— 1=-0.04 mag. The black dotted line represents the 
photometric SSP model predictions computed by adopt- 
ing temperature-gravity-metal licity relations from extensiv e 
photometric stellar libraries (|Vazdekis et al.l 1 19961 . l2010(h 
The colour predictions resulting from these two approaches 
give very similar results (with differences of the order of 
0.01-0.02 mag) showing the reliability of the flux-calibration 
of the model SEDs (see Paper I for more details). 



Figure [T] also shows that the models reproduce the ob- 
served colours fairly well. Note that the V— I colour is par- 
ticularly relevant to test the reliability of our models. This 
colour is heavily sensitive to the details of the method used 
to join the three libraries (MILES, Indo-U.S. and CaT). 
Indeed, the V magnitude is entirely measured within the 
MILES range, whereas the I filter response falls within 
the Indo-U.S. /CaT range. The good match to the observed 
colours shows the robustness of the method and the flux- 
calibration of the MIUSCAT SEDs. It is also remarkable the 
ability of our models to reproduce the B— V colour. This has 
been for a long time a problem for stellar popula tion syn- 
thesi s models. As quo ted in previous works (e.g. IWorthevI 
1 1994 iMarastonl [20051 ). SSP model SEDs tend to provide 
systematically redder B— V colour than the observed values 
(~0.05 mag). One of the reasons of these discrepancies may 
be ascribed to difficulties of the theoretically-based models 
in modeling the spectral energy distribution around ~ 5000- 
6000 A, where the theoretical spectra predic t a flux excess. 
The use of empirical libraries (Pickles 1998) has proven to 
be effective at improving the match with colours of luminous 
red galaxies (|Maraston et alj 2009). Likely, the use of em- 
pirical libraries with good flux-calibration quality and good 
parameter coverage, including the metal-poor regime, has 
allowed us to improve the matching of the GC colours in 
comparison to previous works. The improved match of the 
B— V vs [Z/H] relation when using MILES-based models in- 
stead of the STELIB-K urucz library has been also shown in 
iMaraston et all l|201ll) . 

In Figure we show the results obtained for the U— V 
and U— B colours. We found that colours involving the U 
band are significantly more difficult to reproduce. Indeed, 
our models predict bluer colours («0.1 mag in the high 
metallicity regime). As explained in Sect. 12. ll the U magni- 
tude is not fully derived from the MIUSCAT SEDs, which 
do not extend to the blue end of the U filter respon se. The 
missi ng flux is accounted for with the aid of the (|Pickled 
1998) stellar library. However we rule out this calibration as 
the a main source for the obtained discrepancy, as the ob- 
tained magnitudes following two different approaches (the 
second one employing SSP model SEDs from other authors) 
do agree within ~0.02 mag. Moreover, even the photomet- 
ric predictions show a similar level of disagreement with the 
data. A possible explanation may be linked to the various U 
filter response curves that can be adopted to reproduce the 
observations. Figure [2] shows the obtained results for vari- 
ous U filter responses, which differ from our reference filter 
IIBuser fc Kurucdll978T ). Whereas for other Johnson broad- 
band filters the adoption of different references for a given 
filter usually implies minor changes, in the case of the U 
filter, the exact definition of the curve makes a significant 
impact on the obtained colours. For instance, if we use the 
Johnson (1965) response curve, we obtain a much better 
match. Since it is difficult to trace back the compilation of 
data to recover which filter has been used for every GC, the 
U— V and U— B colours are subject to this source of uncer- 
tainty. 

4.2 M31 globular clusters 

In this section we compare our models with the glob- 
ular clusters of M31. For this purpose we use dif- 
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Figure 1. Optical colours in several bands as a function of metallicity for the MW globular clusters (Harris 1996 catalog; 2010 edition). 
Filled (open) circles represent glob ular clusters with reddening E(B— V)^0.5 (>0.5), de-reddened according to the extinction law of 
ICardelli, Clayton, fc Mathij jl989h . The magenta line shows a second order polynomial fit to the data with E(B— V)^0.5. The black 
dotted lines show the photometric predictions of the MILES models. The synthetic colours derived from the MIUSCAT SSP spectra for 
models with age of 12.6 and lOGyr, for different metallicities and Kroupa Universal IMF, are shown in solid and dashed green lines, 
respectively. 
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Figure 2. The same as figure [T] but for U— V and U— B colours. The green lines show the synthetic colours derived by using different 
definitions of the U filter: solid for the Buser & Kurucz (1978) filter, long-dashed for the Johnson (1965) filter, dotted for Bessel (1990) 
and short-dashed for Fernie (1974). 



ferent photometric catalogs and metallicity estimates. 
We consider t he spe ctroscopic metallicity determinations 
of |Fan et alJ l|2008l). which is a compilation of vari- 



ous, source«_J^^chra^Brodie^&^en3 Il99ll ; iBarmbv et aL 
l200d : iPerrett et all 120021 ). the Revised Bologna Catalo£ 
(RB C), which are det ermined in a homogeneous man- 
ner dGalleti et all [ 2009) and the spectroscopic catalog of 
ICaldwell et all l|201ll ). hereafter Cll. Integrated colours in 
the Johnson-Cousins system have been taken from the RBC, 
while photometry in the SDSS bands is from lPeacock et al.l 
(2010). The extinction in each band has been derived from 
the red dening values given in the corresp onding catalog 
and the ICardelli. Clayton, fc Mathisl (Il989l ) extinction law. 
Given the higher uncertainty in the M31 measurements com- 
pared to the MW GCs, we only consider clusters with low 
extinction values (E(B— V) ^0.5) to avoid unreliable colour 



estimates. Only globular clusters classified as old in the cat- 
alogs are considered. 

Figure [3] shows the Johnson-Cousins colours of the M31 
GCs as a function of metallicity. Given the high dispersion 
of the points we found that the polynomial fit is not a good 
description of the data. To guide the eye, we overplotted the 
median of the colours computed in intervals of metallicity 
and compare it to the fit to the MW GCs. We also show the 
model lines corresponding to the SSPs with ages ranging 
from 6 to 12.6 Gyr (Kroupa Universal IMF). The Vega sys- 
tem used to compute the synthetic magnitudes is the same 
as adopted for the MW data. The SDSS colours in the AB 
system are shown in Figure [3] Note that the scatter ob- 
tained for the M31 data in these two figures is considerably 
larger than that observed in the MW clusters. Moreover, as 
shown by the arrows, the M31 clusters display a considerable 
amount of reddening. 



http: / /www. bo. astro.it/M31 



Although the absolute colours of M31 can be subject 
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Figure 3. Optical Johnson-Cousins colours plotted a gainst the meta llicity for M31 globular clusters. In the first row of panels the 
metallicity and reddening mea surements are taken f rom lFan et al1l2008l , in the second row they come from the RBC, whereas in the last 
row we consider the catalog of lCaldwell et al. (2011). Only clusters with reddening E(B— V)sJ0.5 are shown. The vertical arrow indicates 
the mean amount of de-reddening applied. The magenta thick solid line indicates the median to the data calculated in metallicity bins of 
—2, —1.5, — 1, —0.5, and 0.5 dex. The green lines are the synthetic colours derived from the MIUSCAT SSP SEDs for Kroupa Universal 
IMF and different ages: 12.6 (solid), 10 (dashed), 8 (dotted) and 6Gyr (dashed-dotted). 



to several uncertainties, a robust result inferred from Fig- 
ures 3 is that the slope of the colour-metallicity relation 
is flatter in the M31 clusters. Given that the MW slope is 
consistent with that of the old models, we argue that in 
M31 a significant age variations with metallicity is needed 
to account for the different slope. A flatter relation than 
that indicated from the models is also seen in Figure 4 with 
SDSS colours. Old stellar populations in the lower metallic- 



ity regime and young ages for the metal rich clusters could 
explain such a trend. We obtain the worse fits to the data 
with the old models for the V— R and g— r colours. Indeed, 
the M31 GC colours are in general bluer than their Galac- 
tic counterparts, mainly in the high metallicity regime, as 
it can be inferred from comparing to the MW fit. The rea- 
son for this discrepancy at nearly solar metallicity may be 
linked to the fact that M31 clusters show higher metallicities 
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than their MW counterparts (e.g.. iBeaslev et "al1l2005l ). In 
general the SDSS colours, even for intermediate metallici- 
ties, are matched better with younger models. However, it 
is worth noticing that the choice of the catalog may change 
our conclusion. In fact yo ung models are r equired to fit the 
trend in the g—r colour of lFan et al.l l|2008l ) dataset, whereas 
models of 6-8 Gyr provide a fair match for all the colours of 
the RBC and Cll catalogs. From the comparison with our 
model lines the best match to the data is obtained combin- 
ing old models at lower metallicities and younger models at 
the high metallicity regime. 

In IPeacock et al.l J201l | ') the syntheti c q—r colour de- 
rived from the models of Vazdekis et al.l (|2010l ) has been 
compared to the SDSS colours of M31 clusters in a sim- 
ilar way as it is done here. Since the g—r colour falls in 
the MILES range, it coincides with the predictions of the 
MIUSCA T models presen t ed her e. In IPeacock et al.l (120111 ) 
paper the IVazdekis et al.l l|2010h colours turned out to be 
offset by ~0.1 mag at solar metallicity with respect to M31 
GC colours. Although for the oldest model considered here 
we find similar difference in colour, by using different colours 
and different photometric systems we conclude that the data 
ca n be matched b y con sidering younger models than those 
in IPeacock et al. (2011). Finally, as we obtain with the old 
models very good fits to the MW clusters, which have more 
reliable data, and since the M31 cluster colours in the same 
photometric bands differ from the ones of the MW, we can 
conclude that we require younger models, at least in the high 
metallicity regime, to fit the data. 

Our results agree with evidences supporting that metal- 
rich GCs tend to be younger than t he more metal-poor 
clust er s, at least in t h e MW /e.g.. | Marfn-Franch et all 
20091; |Pe Angeli et all 120051 ; iBarmbv fc Huchral |200oF 



Rosenberg et al.lll999h . Moreover, it has been claimed that 



the GCs of M31 tend to be on average yo unger than their 
Galactic counterparts l|Beaslev et al.ll2005l ). though this re- 
sult is yet to be conclusive. Spectroscopic ages for the 
high-metallicity M31 clusters have been derived in Cll. In 
the low-metallicity regime the authors could not get ac- 
curate age determination, as these clusters are contami- 
nated by blue HB stars, hence we can not test the age 
variation with metallicity. For the high-metallicity clusters 
in our sample ([Z/H]> —1) the Cll determinations give 
a mean age of 11.3 Gyr, that apparently contradicts the 
young nature of the M31 clusters. Note however that these 
authors employed the Balmer absorption line indices as 
age indicators and it is well known that this method can 
give very uncert ain absolute ages, even larger than the age 
of the Universe (I Vazdekis et all 1 200 ll; ISchiavon et al.l l 2002j ; 
iMendel. Proctor fc Forbes! 120071: iPoole et al.ll2010l ). More- 
over, the models of lSchiavon 1 2007 ) used in Cll give broad- 
band colours in excellent agreement with our values. For in- 
stance, a solar metallicity model with age of 11.2 Gyr has a 
B— V colours of 0.978 mag, very similar to B— V=0.977 mag 
that we obtain for the MILES models. Both photometric pre- 
dictions are much larger than the M31 colours (Figure [3}, 
confirming the discrepancy between the spectroscopic and 
photometric age measurements. 

Furthermore, given the high uncertainties involved in 
the extinction, even our results should be taken with cau- 
tion. Indeed, the parameter that varies mostly among the 
three catalogs is the reddening. The Cll estimates tend to 
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Figure 5. Redshift distribution for the LRG sample, which is 
composed of 168737 galaxies. 



be biased towards larger values than the other two catalogs, 
leading to bluer de-reddened clusters (see for instance the 
g—r colour in Figure[4]|. In Cll the quoted uncertainties for 
the E(B — V) is 0.1 mag, hence of the order of colour vari- 
ations due to different ages. Therefore, it is not possible to 
draw robust conclusions. 



5 COLOUR EVOLUTION OF LUMINOUS RED 
GALAXIES 

In this section we consider the Lumin ous Red Galaxy 
(LRG) sample of lEisenstein et all (|200ll ). for which typi- 
cally the definition of early-type galaxy (ETG) also applies. 
The LRG sample has been extracted from the SDSS spec- 
trosc opic galaxy survey (Data Release 7, lAbazaiian et al.l 
l2009h . by selecting galaxies with the GALAXY_RED flag 
in the SDSS database. LRGs constitute a uniform sample 
of objects with the reddest rest-frame colours, selected with 
cuts in the (g—r , r—i, r) colour-colour-magnitude cube (see 
lEisenstein et alj|200ll ). To obtain the same galaxy popu- 
lation at different redshifts, the LRG sample has been se- 
lected using two different selection criteria below and above 
z=0.4, as the Balmer break moves from the g band to the r 
band at this redshift. To increase the redshift range, galax- 
ies from the 2df SP SS LRG and Quasar (2SLAQ) survey 
ijCannon et al]|2006l ) have been also included (~15000 ob- 
jects) and they contribute to the high redshift tail of the red- 
shift distribution (see Figure[5]). By restricting the sample to 
galaxies with reliable redshift determination (zWarning=0) 
we end up with ~170000 galaxies. The SDSS database pro- 
vides a variety of measured magnitudes for each detected 
object. Throughout this paper, we use dereddened model 
magnitudes because they provide an unbi ased colour esti- 
mate in the absence of colour gradients (jStoughton et al] 
l2002h . 

Previous works l)Eisenstein et al ] l200ll : IWake et al.l 

2006) have highlighted the difficulties of stellar popula- 
tion models in reproducing the colour evolution in observed 
frame of LRG at intermediate redshift (0.1 < 2 < 0.7), with 
the prediction of too red colours in the g—r and too blue 
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Figure 4. Same as in Figure \3\ but for colours in the SDSS bands. The errorbars in the lower-right corner indicate typical errors on the 
metallicity and the photometric measurements. 



in the r—i for the lowest redshift LRGs, while the opposite 
effect is observed for the highest redshift LRGs. Similar dis- 
agreements have been reported for red sequence colo urs in 
galaxy clusters at similar redshifts (|Wake et al.l l'2005'). 

In Figure[6]we compare the predictions from the MIUS- 
CAT models to the colour evolution of the LRG sample. Ac- 
cording to the fiducial model for the formation of LRG, in 
the simplest case we assume that they formed in a uniform 
burst at high redshift (z=5) and evolve passively thereon. 
Solar metallicity and Kroupa Universal IMF are assumed. 
The observed colours are reproduced by shifting the model 
spectra to the corresponding redshift and taking into ac- 



count the ageing effect. Synthetic colours in the SDSS fil- 
ters have been computed in the AB system (|Oke &: Gunnl 
1983). For consistency, the observed colours from the SDSS 
database have been transformed to AB system by apply- 
ing the relative shifts: uab = usdss — 0.04 mag; zab = 
zsdss + 0.02 ma£|. The model predictions are shown for 
the redshift interval where the MIUSCAT spectra cover the 
r— i and i— z colours. The g filter falls outside the measur- 
able range at z=0.05, thus colours containing u and g are not 
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Figure 6. Colour evolution of the LRG sample in the r—i (left-hand panel) and i— z (right-hand panel) bands. The contours enclose 
20%, 40%, 60% and 95% of the galaxies. The red thin line shows the median colour at different redshift bins equally spaced in intervals of 
0.025. The black thick line shows a passive evolution model with solar metallicity and redshift of formation of 5 (corresponding to an age 
of 12.6 Gyr at z=0) for a Kroupa Universal IMF. The dot-dashed line indicates a model that includes a 5% contribution in mass from a 
metal-poor ([Z/H]=-2.32) stellar population. The residuals between the synthetic and observed colours are shown in the bottom panels 
using the same scale. The bottom panels show the colour evolution for a subsample of LRGs having absolute magnit ude: Mj < -23 mag. 
The do tted vertical lines indicate the redshift range where these bright galaxies are a coeval population according to iToieiro fc Percivall 
(2011bl) . 



shown. We see that the simplest model is able to reproduce 
the observed distribution at low redshift in both colours. The 
model starts to deviate from the median colours at z~0.4, 
producing too red colours. 

In Figure [6] it is represented also a second model, 
where a small fraction of the mass (5%) is contributed 
by a metal-poor population ([Z /H]=-2.32). It is expected 
|Worthev et alj|2005l : iPagelll 19971 ) that galaxies exhibit some 
spread in the metallicity distribution with a contribution 
also from metal-poor stars. The effect of including stars with 
such a low metallicity is to bluen the colours, with the result 
of loosing the match with the observed distribution in the 
lowest redshift range, but improving it at high redshift. It is 
worth to say that the differences between this model and the 



median of observed colour is only w 0.02 mag. The uncer- 
tainties derived from the zero-point calibration of the SDSS 
system (0.01 mag) and from the choice of the IMF (0.01 
mag of difference between Salpeter and Kroupa Universal) 
could also account for this difference. Hence, it turns out 
to be difficult to discriminate among these two possibilities. 
Nevertheless, in both cases the MI USCAT models perform 
rather well, confirming the results of lMaraston et all (2009). 
according to which models constructed with empirical stellar 
libraries provide a good fit to the LRG colour distribution. 

Although the LRG sample has been selected to follow 
a passive evolution model (see lEisenstein et alJ feOQl). this 
does not guarantee that all the galax ies accomplishing the 
selection criteria follow this model. In lToieiro et alj (|2011al ) 
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redshift 

Figu re 7. Redshift distribution for the ETG sample of 
iNair fc Abraham! teulCh . The red lines indicates our sample at 
z~0.04, which includes 370 objects. 

it has been shown that a non negligible fraction of LRGs 
does not evolve passively in a dynamical sense. They can 
undergo mergers and the sample can be contaminated with 
galaxies that belong to the LRG sample for a short period 
of time. However, the brightest tail of the sample turns out 
to be coeval (i.e. with almost no merging) for galaxies in the 
redshift range: 0.22-0.45. For this reason, in the bottom pan- 
els of Figure [6] we compare our passive evolution models with 
these bright galaxies (Mi < -23 mag). For this subsample, 
the high-redshift galaxies are truly the progenitors of the 
low-redshift ones. In the redshift range where the galaxies 
appear coeval, our passive model performs remarkably well. 



6 NEARBY EARLY-TYPE GALAXIES 

A further insight into the capabilities of the models can be 
obtained by comparing them to a sample of nearby galaxies, 
for which we can use three different colours as constraints. 
For this purpose we select a homogeneous sample in the 
Local Universe with a large numb er of SPSS galaxies wit h 
visual morphological classification l|Nair fc Abrahamll2010l ). 
We select a subsample of early-type galaxies (morpholog- 
ical type T ^ —4) with reliable redshift measurements. 
We further restrict our sample to a narrow redshift slice, 
0.035 < z < 0.045, to end up with 370 targets. This choice 
of redshift is motivated by two main reasons. First, at this 
redshift the number of galaxies peaks (see Figure [7J. Sec- 
ond, given the limited spectral coverage of the MIUSCAT 
models in the blue end, the correction factors we need to 
apply to measure the SDSS u magnitude become too large 
for redshifts above this value, as the wavelengths at which 
the u filter response peaks start to fall outside the spectral 
range of the models (see Paper I). Before comparing the 
colours derived from the MIUSCAT SEDs, once redshifted 
to z=0.04, we have tested that the selected redshift interval 
is small enough that colour differences due only to variation 
in redshift within such interval are negligible (<0.01 mag). 

As for the photometric estimates, we rely on the model 
magnitudes from the DR7 catalog. It is worth noticing that 



all the colours involved in this study belong to the optical 
spectral range, which are known to be heavily affected by the 
age/metallicity degeneracy, irrespective of the colour-colour 
diagnostic diagram in use. Therefore the results obtained in 
this section are not intended to provide a well-constrained 
picture for the stellar populations of these galaxies, which 
might require a deeper analysis involving colours in other 
wavelengths and/or spectra. 

In Figure |8] we show the colour distribution of the 
ETG sample compared to the MIUSCAT SSP models with 
ages ranging from 3 to 14 Gyr for three metallicity values: 
[Z/H]=-0.4, and 0.22. We also include an arrow indicating 
the effect of internal dust extinction on the observed colours. 
We assume Av = 0.2 mag, in agree ment with the typica l 
value derived for low redshift LRG (|Toieiro et al.l l2011al ) . 
The de-reddening in the various SDSS bands has been 
computed with the e xtinction law: t\ — tv(A/5500A) -0 ' 7 
|Charlot fc Fal]|200Ch . We see that none of the SSP models is 
able to match the observed colour distribution. Although the 
fiducial model for the formation of ETGs predicts that they 
form in a uniform burst at high redshift with high metallic- 
ity, the old models tend to produce too red colours in the 
g— t vs u—g plot. On the contrary, the r— i colours turn out 
to be too blue with respect to the real galaxies. Hence, this 
simple picture does not seems to appropriate to describe the 
stellar populations of this galaxy sample. 

The MIUSCAT models adopt the stellar isochrones of 
iGirardi et ail (|2000h (i.e. Padova 2000), as described in 
Paper I. It has been claimed that t hese tracks provide 
blue c olours for old elliptical galaxies. iBruzual fc Charlotl 
(2003 ) claimed that m odels computed with the isochrones 
of iBertelli et alj <| 19941 ) (i.e. Padova 1994) are preferred be- 
cause redder V-K colours are predicted, in better agreement 
with the observations. For the optical spectral range the net 
effect of adopting these two sets of stel lar tracks is almost 
neglig ible. Using the Padova 1994 tracks. [Bruzual fc Charlotl 
(2003) obtained redder colours than using Padova 2000: u—g 
by ~0.02 mag and g—r and r—i by ~0.01 mag. Therefore 
these results tend to worsen the obtained mismatch for the 
bluer colours. 

To understand better the origin of this mismatch and 
to avoid possible sample selection biases we have divided 
the ETG sample in intervals according to their masses. In 
Figure [9] we show the colour distributions corresponding 
to galaxies with velocity dispersions above and below 180 
km s _1 . This value allows us to separate galaxies with dif- 
ferent properties, while keeping almost equally populated 
these two galaxy mass ranges. ETGs with high velocity dis- 
persion are generally found to exhibit older ages and higher 
metallicities than their low velocity dispersion co unterparts 
(|Trager et al.ll200d : fSanchez-Blazquez et alj |2006). resulting 
in redder colours. We see that none of the SSP models plot- 
ted in this Figure is able to match the distribution of massive 
ETGs. Galaxies with lower velocity dispersion tend to be 
matched with models of lower metallicity and age ~8 Gyr, 
even though the predicted u—g remain too blue (~0.1 mag). 
Although a better agreement is obtained, the models do not 
provide fully satisfactory fits. 

The colour predictions computed by adopting SSP 
models from different authors are shown in F igure 1101 
We considered the BC03 and iMarastonl (|2005l ) models, 
the two based on the BaSeL theoretical stellar library 
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Figure 8. Observed colour distributions in the SDSS bands for the ETG sample at z~0.04. The contours enclose 20%, 40%, 60% and 
95% of the galaxies and the cross indicates the median colours (u — g = 1.83, g — r = 0.83, r — i = 0.40 mag). The arrow indicates the 
colour de-reddening when an internal extinction of Av = 0.2 mag is assumed. Synthetic colours derived from the MIUSCAT SEDs for a 
Kroupa Universal IMF are overplotted. The different lines join models with the same metallicities for ages ranging from 3 to 14Gyr. 



Leieune. Cuisinier, fc Buser 19981') and t he [Maraston et al.l 
2009) mode ls that use the | Picklesl (Il998l) empirical li brary . 



The BC03, iMarastonl (|2005l ) and iMaraston et all 
adopt a Salpeter IMF, whereas the MIUSCAT models are 
shown for both Salpeter and Kroupa IMF, showing that the 
choice among these commonly used IMFs has little effect on 
these colours. Although none of the models can fairly re- 
produce the observed distribution, it is remarkable the ad- 
vantage of using models based on empirical-based libraries. 
These models predict g— r bluer by 0.06 mag and r—i redder 
by 0.05 mag compared to the mod els based on theoretical 
libraries. The lMaraston et alj (|2009T l models predict slightly 
bluer g—r and redder r—i colours than the MIUSCAT SEDs 
by «0.01 mag, but also redder u— g colour. 

6.1 IMF effects 

Here we explore how these colours vary when adopting dif- 
ferent IMF shapes and slopes. We consider the power-law 
unimodal IMF defined in Vazdekis et al. (1996) for different 
slopes for the solar metallicity SSP models with the same 
ages considered in Figure [8] Figure [11] shows a very modest 
variation of the bluer bands. However varying the IMF turns 
out to have a significant effect on the r—i colour. Indeed, by 
assuming a very steep IMF slope (/i = 2.3) the colour distri- 
bution in the g—r vs r—i diagram can be matched, though 
for relatively young ages (3-4 Gyr). On the other hand the 
g—r vs it— g colour distribution is not reproduced for none 
of the represented models. 



with a primordial chemical composition ([Z/H]=-2.32) and 
very old age (14 Gyr). This should resemble the first popula- 
tion of stars formed before being enriched by the subseq uent 
star formation. As shown also bv lMaraston et all (|2009l ). the 
contribution of metal-poor stars is effective in bluening the 
g—r, without perturbing too much the r—i colour. The re- 
sulting MIUSCAT two-burst colours are shown in Figure 
1121 As expected, all the colours become bluer, improving 
the match with the observed colour-colour distribution in 
the g—r vs u—g diagram. The resulting r—i colour becomes 
slightly bluer, but since the effect of the metal-poor compo- 
nent is stronger on the bluer bands the colour distribution 
in g— r vs r—i improve s as w ell, in agreement with the find- 
ing of lMaraston et al] (|2009T ) . However the inclusion of such 
metal-poor component alone is not sufficient to fairly repro- 
duce the observed colour-colour distributions. 



6.3 Two-burst model: inclusion of a young stellar 
population 

A similar effect is seen when a young component is included. 
In Figure[l3]the plotted model includes a small contribution 
(3% of the total mass) of a young stellar population of 1 Gyr. 
In this case the effects on all the colours is more pronounced 
than in Section \6. 21 This model can successfully explain the 
colour-colour distribution in the g—r vs u—g diagram, but 
not the r—i colour, since the predicted values are far too 
blue («0.06 mag). Therefore, these models are not able to 
match the two colour-colour diagrams simultaneously. 



6.2 Two-burst model: inclusion of a metal-poor 
component 

We relax here the assumption that ETGs are SSPs and use 
stellar populations with two main bursts to model these 
galaxies. The first model tested is the same used for the 
description of the LRG colour evolution shown in Section 
5. We include a small contribution of 5% to the total mass, 



6.4 Three-burst model: chemo-evolutionary 
model with varying IMF slope 

Here we switch for the first time to the combination of more 
than two stellar populations. In Sections 16.21 and 16.31 we 
have seen that the combination of an old metal-rich stel- 
lar population with a small fraction of either a metal-poor 
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Figure 9. Upper panels: observed colour distribution for ETGs with velocity dispersion below 180 km s _1 , compared to the same SSP 
models as in Figure [8] Bottom panels: colour distribution for ETGs with velocity dispersion above 180 km s —1 . The median colors of 
the observed distribution for low and high velocity dispersion objects are, respectively: u — g = 1.76, g — r = 0.81, r — i = 0.38 mag and 
u — g = 1.87, g — r = 0.85, r — i = 0.40 mag. 



or a young component has significant effects on the re- 
sulting colour-colour plots. Furthermore in Section 16.11 we 
tested the effects on the colours from non standar d IMFs. 
In fact, very recently Ivan Dokkum fc Conrovl | |20l4 showed 
spectroscopic evidences pointing towards steeper IMFs in 
massive galaxies, in agreement with previous studies (e.g. , 
Faber fc French! 1 198CI ; ICarter. Visvanathan fc Pickles! 1 19861 ; 
Vazdekis et all 1 19961 ; ICenarro et al.ll2003l ). Therefore it is 
worth to test the colours that result from combining all 
these stellar popu lations. An example of such scenario is 
the one derived in lVazdekis et all (|l997l ) from the best fits 
to both broad-band colours (from U to K) and the Lick line- 
strength indices, of a number of prototype massive galaxies, 
employing a full chem o-evolutionary version of our models 
i|Vazdekis et al.lll996h . 

In this scenario the IMF varies with time, being skewed 
towards high-mass stars during an initial short period 
of time (<0.5Gyr) that quickly enriches the interstellar 
medium and, afterwards, the IMF becomes even steeper 
than the Salpeter one. This is reproduced with a multiple 



bursts model by adding on top of the main episode of star 
formation, that forms 87% of the mass, two populations. In 
the first population 10 % of the mass formed at 14 Gyr with 
a top-heavy IMF (slope ^i=0.8). Half of the mass formed 
in this burst is associated with a metal-poor component, 
[Z/H]=-2.32 dex (as in Section |6.2[) . and the other half has 
subsolar composition ([Z/H]=-0.4). Another 3% of the mass 
is formed at the age of 1 Gyr (as in l6.3|) . The stars formed in 
the main episode and at late times are skewed towards low- 
mass (/x=2.3) (as in Section ^, l[l . A bimodal IMF is assumed 
for all these models. The adoption of the IMF with slope 
2.3 in the main burst is consisten t with recent evidences 
for a steep IMF in massive galaxies |van Dokkum fc Conrovl 
|2010D . 

This model is shown in Fig. 14. For solar metallicity it 
is quite successful at matching the bulk of the galaxy distri- 
bution, although some problems persist in the r— 1 colour, 
which remains too blue. Although quite unlikely, part of this 
discrepancy can also be attributed to internal extinction, 
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Figure 10. Same as in figure[8]for stellar populations models from different authors. For this comparison only models for solar metallicity 
are shown. The blue line corresponds to the BC03 models assuming a Salpeter IMF. Red dashed line represents the M05 models based 
on theoretical stellar atmospheres, w hereas the solid red line shows the models from the same authors based on the Pickles empirical 
stellar library llMaraston et al] [2009). In these two cases a Salpeter IMF is adopted by these authors. Finally the black lines indicate the 
MIUSCAT model tracks adopting the Salpeter (solid line) and Kroupa Universal (dashed line) IMF. 




Figure 11. Same as in Figure \E\ for unimodal IMF with different slopes. 



leading too red colours that we cannot fit with dust-free 
models. 



6.5 Three-burst models: Montecarlo simulations 

In order to explore all the possible combinations of SFHs 
that can match the two colour-colour distributions we have 
run a set of Montecarlo simulations, where we allow all 
the stellar population parameters to vary: age, metallicity 
and IMF slope. We only consider seven representative IMF 
slopes: 0.3, 0.8, 1.3, 1.8, 2.3, 2.8 and 3.3. To keep the case as 
general as possible we use three bursts, where the fraction 
of mass involved in each of the burst is also a free parame- 
ter. Figure [Pol shows the results for 100000 simulations. On 
one hand we see that the cloud of models does not properly 
match the observed galaxy distribution in the g—r vs u— g 



plot. On the other hand, the cloud of models matches the 
region enclosed by data in the g—r vs r—i plot. 

To understand which kind of SSP mixtures most likely 
resemble the data we selected the closest 1700 points to the 
median of the observed colours in each of the distributions. 
Among these combinations we only choose those that are 
common to these two distributions to end up with 26 models, 
whose SFHs are shown in the Appendix. We have ordered all 
these models following an increasing contribution of young 
stellar populations, i.e., with ages smaller than 5 Gyr. We 
see that a large number of models require very important 
contributions from SSPs with old ages, high metallicity (so- 
lar or higher) and IMF slopes steeper than Salpeter (i.e. 
ft > 1.3). All these models require smaller, but not negligi- 
ble, contributions from young or/and metal-poor stellar pop- 
ulations. Finally SSPs with slopes lower than the Salpeter 
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Figure 12. Same as in Figure [8] but for a two-burst model that includes a metal-poor population ([Z/H]=-2.32) of 14Gyr, which 
contributes with a 5% to the total mass. 




Figure 13. Same as in Figure 1121 for a two-burst model including a 3% contribution in mass from a 1 Gyr stellar population. 



value seem to be required for a number of models. A mean 
solution for these models can therefore be summarized as 
mainly composed of metal-rich and dwarf-dominated stellar 
populations with much smaller contributions from young or 
metal-poor components. Interestingly the scenario discussed 
in Section T6.4I would be a good representative of this mean 
trend of solutions. In fact we nearly recover it in the models 
shown in the first and second rows of panels of Figure IA1I 
Whereas many studies support the presence of old ages and 
metal-rich stellar po pulations in m assive galaxies (see for ex- 
ample the review of Renzini 2006), the requirement of non 
st andard IMFs is a matter of deb ate as recently pointed out 
bv lvan Dokkum fc Conrovl (|2010h . Among our solutions, the 
steeper IMF slopes are recovered for the highest metallicity 
objects. In Figure [161 we show the mass weighted mean IMF 
slope as a function of the mass weighted mean metallicity for 
the 26 solutions. The existence of a I MF slope - me t allicit y 
relation has already pointed out by ICenarro et all (2003). 
Interestingly, the linear fit to the simulated points lies close 



to the relation derived by the above authors by means of 
near-IR IMF sensitive spectral features. It is worth to note 
that, using completely independent diagnostic techniques, it 
is inferred a similar IMF-metallicity relationship in the sense 
that larger metallicities correspond to larger effective IMF 
slopes. 

There are a number of models that depart significantly 
from the mean trend depicted above (see last rows of Fig- 
ure lAl[) . These models require important contributions from 
stellar populations with ages smaller than ~5 Gyr. However 
we also see that most of these models also require steeper 
IMF slopes than the models that are clearly dominated by 
old stellar populations. This can be explained by the fact 
that the colours redden with increasing age, metallicity or 
IMF slope as shown in Paper I. Such degeneracy between 
these three parameters for the optical colours does not allow 
us to provide a unique solution. To constrain the solutions it 
is required to use spe ctroscopic data or/an d colours in other 
spectral ranges (e.g., IVazdekis et"aflll997l ). 
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Figure 14. Same as in Figure[8]but for the SFH derived from the chemo-evolutionary model described in Vazdekis et al (1996, 1997). 
We model it by adding on the top of the main episode of star formation the contributions of two stellar populations: 10 % of the mass 
formed at 14Gyr with a top-heavy IMF (slope fi=0.8). Half of the mass formed in this burst is in metal-poor component ([Z/H]=-2.32 
dex) and half of it has subsolar composition, [Z/H]=-0.4 dex. Finally, there is a 3% contribution to the total mass from a component of 
1 Gyr. The stars formed in the main epis ode and afterwa r ds are skewed towards lower masses according to a steeper IMF slope (fi=2.3). 
We use here the Bimodal IMF defined in Vazd ekis et al.l 1 1199611 . 




Figure 15. The observed colour-colour distributions are compared with simulated points derived from the three-bursts Montecarlo SFHs 
described in Section 16.51 



However it is worth recalling that none of the SFHs ex- 
plored in the simulations provide satisfactory fits to these 
two observed colour-colour distributions. Indeed, the depar- 
ture of the models from the median of the observed colours 
is non negligible, with the u — g colour showing the largest 
residuals. It is beyond the scope of this paper to recover the 
SFHs of the SDSS early-type galaxy sample, but the fact 
that in these complex models we allow to vary all, the age, 
metallicity and IMF slope of the SSPs, without reaching a 
fully satisfactory solution, supports the need for a fourth pa- 
rameter to be taken into account in the models. We consider 
such a possibility in the next section. 



6.6 Effect of a-enhancement 

As explained in Paper I the MIUSCAT SSP SEDs employed 
here are computed on the basis of scaled-solar isochrones and 
stellar spectra with nearly scaled-solar abundance ratios for 
the high metallicity regime that is required to fit massive 
galaxies. However, it is well known that the spectra of these 
galaxies show an enhancement in [Mg/Fe]. To assess the ef- 
fects of Q-enhanced element partitions on the SSP colours 
we correc t them, in a differe ntial way, with the aid of the 
models of ICoelho et alj ((2007) . These fully theoretical mod- 
els are constructed o n the basis of a hig h resolution library 
of synthetic spectra l|Coelho et al. 11200a ). where all the a el- 
ements, including magnesium, are enhanced with respect to 
iron with the same proportion, i.e., [a/Fe]=0.4. 
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Figure 17. The observed colour-colour distributions are compared to the MIUSCAT SSP models once corrected from the a-enhancement 
effect. There are three shaded regions corresponding to each metallicity. For each of these regions the upper envelope indicates the 
MIUSCAT SSP models (same as in Figure [8), where as the lower e nvelope shows the resulting models once corrected from the a- 
enhancement effects with the aid of the models of lCoelho et al. | J2007T) . For this purpose we estimate the colour difference obtained from 
comparing the scaled-solar and a-enhanced ([a/Fe]=0.4 dex) models from these authors. 
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Figure 16. Mean IMF slope versus mean metallicity for the 
best solution SFHs from the Montecarlo simulations described 
in Section 16.51 The solid black line shows the relation found in 
ICenarro et al.l (2003), whereas the dotted line is the linear fit to 
the points, obtained with a least squares fit in the X direction. 



We correct our MIUSCAT SSP colours with the dif- 
ference in colour obt ained from th e scal ed-solar and a- 
enhanced models of ICoelho et al.l (|2007l ). keeping con- 
stant the iron abundance. The corrected SSP colours are 
shown for different metallicities in Figure 1171 Interest- 
ingly, the effect of a-enhancement tends to bluen the u—g 
and g—r colours and to redden r— i, as required by the 
data. Note that although this effect is relatively modest 



it certainly goes in the right direction for improving our 
fits to the most massive galaxies of our sample. Prelim- 
inary tests performed with our own version of the mod- 
els constructed with a-enhanced mixtures based on the 
M ILES library, making use of the [Mg/Fe] de terminations 
of lMilone. Sansom. fc Sanchez-Blazaued l)201lf ) , lead to even 
larger differences with respect to the scaled-sol ar predictions 
in the same directions as with the models of ICoelho et al.l 
(|2007h . 



7 DISCUSSION AND CONCLUSIONS 

In the present work we have explored the constraints on the 
MIUSCAT models provided by photometric data of globular 
clusters and quiescent galaxies. In Paper I we have presented 
the new models and described the improvements. Here, we 
have focused on the comparison of the model predictions 
with broad-band photometry. 

The integrated colours of MW globular clusters can be 
remarkably reproduced by the MIUSCAT models in all the 
visible bands. Only colours involving the U filter can not be 
used to constrain the models, as the uncertainty relative to 
the filter response used to reproduce the U magnitude is too 
high. The ability of our models for matching the colours of 
the MW stellar clusters in the other bands is remarkable. 
In particular, the B— V colour has always challenged models 
based on theoretical stellar libraries, which fail at matching 
the range covered by the V filter. Here we have shown that 
the problem can be solved by the use of empirical stellar 
libraries (see also iMaraston et al.ll201lh . 

The comparison with M31 globular clusters is less 
straightforward but it allows us to test the models even 
in the region of solar metallicity. The comparison with our 
models suggests that the M31 globular clusters are on av- 
erage younger than usually assumed for the MW globular 
clusters. Although an old model (12.6 Gyr) can fit the metal- 



16 Ricciardelli et al. 



poor clusters, we need ages in the range 6-10 Gyr to match 
the colours of the clusters with high metallicities. 

The spectral coverage of the MIUSCAT library allows 
us to test here the models in several bands and, by relaxing 
the assumption that all the M31 GC are old, we can obtain 
a good fit to the data if we assume that metal-rich M31 clus- 
ters are younger than metal-poor ones. The fact that in MW 
clusters we do not need young models may be ascribed to 
the lack of metal-rich clusters. Therefore o ur results disagree 
with previous works (jPeacock et al.|[201lh claiming that the 
g— r colours of M31 globular clusters could not be matched 
accurately with our models. Nevertheless, the large amount 
of extinction in these clusters makes their colours very un- 
certain and not suitable for calibration purposes. 

When the models are compared to the colour evolution 
of LRGs, under the assumption of a passive evolution model, 
the observed colours are fairly reproduced. Both, a simple 
model where stars formed in a single burst at high redshift, 
and a model that includes a small component of metal-poor 
stars, can be used to fit the observed distribution with in- 
creasing redshift. However, in the high redshift interval (z > 
0.4) where the observed colours probe the rest-frame g— r, 
the models predict redder colours than the median distribu- 
tion. 

Although the colour evolution of the LRGs can be fairly 
described within the fiducial model of purely passive stellar 
evolution, the combination of three colours, u— g, g—r, r—i, 
poses a challenge. We find that the colour-colour diagrams 
based on these colours allow us to discriminate among differ- 
ent models. Given the wavelength coverage of the MIUSCAT 
models, we restrict our study to a limited redshift range, at 
z~0.04, where we can measure all these three colours. Our 
study shows that single-burst models do not provide a good 
description for the local population of quiescent galaxies. Ir- 
respective of their metallicity, the old SSPs provide too red 
u — g and g—r colours and too blue r—i colour. We also find 
that the SSP model predictions from different authors do not 
match these colour-colou r distributions. S u ch dis agreement 
has also been shown by IConrov fe Gunnl |2010T l. We note 
that the use of only two colours, u—g and g—r does not allow 
to highlight the problem in the r—i colour, which goes in the 
opposite direction. Nonetheless, an important result drawn 
from such a comparison is t hat models based on e mpirical 
stellar libraries (MIUSCAT; iMaraston et al.l l2009h predict 
colours much closer to the observ ed distribution than mod- 
els based on th eoretical spectra l|Bruzual fc Chariot! 120031 ; 
lMaraston|[2005l ). 

We rule out that the inability of the MIUSCAT SSP 
models for providing fully satisfactory fits to the colours 
of LRGs is due to flux-calibration issues. As shown in Pa- 
per I, the MIUSCAT spectra provide colours with photo- 
metric uncertainties < 0.02 mag. Furthermore, the colours 
of our models are in good agreement with globular cluster 
data. 

We have tested the hypothesis that quiescent galaxies 
formed their stars in a strong burst at high redshift. Our 
study based on the u—g vs. g—r and r—i vs. g—r colour- 
colour diagrams suggest that the constituent stellar pop- 
ulations of local ETGs are not necessarily fully coeval and 
old, and might require small contributions from either young 
or/and metal-poor stellar populations. We have explored 
several models with such complex stellar populations. We 



find that these smaller contributions are very effective at 
bluening the colours and matching the distribution in the 
u—g and g—r colours. Note however that the resulting 1 — i 
colour is too blue. 



The chemo-evolutionary model shown in lVazdekis et all 

(|l996h provides one of the best fits to these two colour- 
colour diagrams. Apart from the inclusion of a metal-poor 
population at early-times and a small burst at 1 Gyr, this 
model adopts a variable IMF scenario. In this model the first 
generation of stars is skewed toward high-mass (top-heavy 
IMF) , while at later times stars formed with a steeper IMF 
slope. The adoption of such a steep IMF is the main rea- 
son for the improvement also in the r—i colour. Although 
this model is just one among several possible solutions, 
it clearly shows the need for a more complex SFH than 
the single burst model to explain the formation of ETGs. 
Moreover, the finding that quiescent galaxies have experi- 
enced a small amount of recent star formation is not new 
and has been revealed by several observational evidences 
dTrager et al.ll2000l; iBressan et al.ll2006l ; iKavirai et alj|2007l ; 
ISarzi et al.ll2008l ; iToieiro et al.ll2011al ) and predicted by hi- 
erarchical models of galaxy forma tion (|De Lucia et al,ll2006l ; 
iRicciardelli fc Franceschinill2010i ). However, it is worth not- 
ing that although the colours obtained from these complex 
models are approaching the observed values, the agreement 
is not satisfactory. 

An experiment done by simulating three-burst Monte- 
carlo SFHs confirmed the impossibility of the three model 
parameters (age, metallicity and IMF slope) to fairly match 
the color plane g—r vs u—g. We are aware of the degen- 
eracies affecting the optical colours, which do not allow us 
to properly constrain the solutions. It is beyond the scope 
of this work to fully constrain the SFHs, for which spec- 
troscopic data or colors in other spectral range would be 
required, but the inability of these simulations to reproduce 
the observed colours suggests the need of another parame- 
ter to model the observations. We have identified that the 
inclusion of a-enhancement improves the fits. A limitation 
of the models based on empirical libraries is that they rely 
on stars from the solar neighborhood, which, most likely, 
has experienced a chemical enrichment that differs from 
that of massive elliptical galaxies. Indeed, it is now well 
established that massive early-type galaxies are enhanced 
in Magnesium with respect to Iron, showing \Mg/Fe] ~ 
0.2 dWorthev Faber. fc Gonzalez! Il992l ; iTrager et all |2000| ; 
lYamada et al. I l2006h . We have estimated the effect of a- 
enhanced models on these colours by applying differentia l 
corrections derived from the models of lCoelho et alj (I2007T ) 
predictions on the MIUSCAT colours. The overall effect is 
in the direction of improving the fits to the observations, as 
the colours in the bluer bands become bluer whereas the r— i 
tends to redden. Therefore, the combination of both effects, 
a more complex SFH and the inclusion in the models of 
abundance ratios in agreement with the observations should 
provide better fits. 

We conclude that the impact of a-enhancement on the 
colours of massive galaxies can not be neglected to properly 
reproduce the observations. Our forthcoming models taking 
into account the a-enhanced mixture of MILES stars will 
clarify this issue. 
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the spectra that are no t observed in the lo cal population of 
massive ellipticals (e.g. iTrager et al.ll2000h . 

Interestingly, within the first group of models we can 
identify a number of solutions that resemble the ones found 
in section l6~4l A small burst at early times occurring with a 
flat IMF and metal-poor stars is followed by a strong burst 
at nearly solar metallicity and steeper IMF and then by a 
small burst at later times (see first and second rows). 



APPENDIX A: RESULTS FROM 
MONTECARLO SIMULATIONS 

In this section we show the SFHs of the best solutions found 
in the Montecarlo simulations described in Section [6.51 We 
selected the 1700 closest points to the median of the ob- 
served colours in each of the colour-colour plots. Then we 
only considered the solutions that satisfy both distributions, 
ending up with 26 models, whose SFHs are shown in Fig- 
ure IA1I The SFHs are ordered according to an increasing 
fraction of mass in bursts younger than 5 Gyr (the first pan- 
els show the solutions with the smaller fractions). 

The resulting SFHs can be divided into two broad cat- 
egories. In the first group we can identify SFHs dominated 
by an old, metal-rich population, enriched with dwarf stars 
(high IMF slopes). On top of the old populations, we find 
small contributions from young populations. In some cases 
we do not see any significant young component. In these 
cases, blue g—r and u—g colours result from the contribu- 
tion of a metal-poor population (see also section |6~2]) . 

The other group of solutions includes SFHs with an im- 
portant component of young populations (^ 5 Gyr). In some 
cases, the fraction of mass formed in the youngest burst can 
be as large as 80% (see SFH in the last row). As already 
shown in Section T6. 31 the age is the dominant parameter at 
decreasing both the g~r and u—g colours and, thus, it pro- 
vides a good match of the g—r vs u—g plot. The SFHs with 
the largest component of young stellar populations also show 
the steepest IMF slopes, which are needed in order to redden 
the r—i colour. This confirms the results found in sections 
6.1, 6.2 and 6.3, where we have illustrated how the IMF is 
mainly affecting the r—i while age and metallicity variations 
are needed to match the g—r vs u—g colour distribution. 

Given the age/metallicity/IMF degeneracies affecting 
the optical colours, the two kind of SFHs produce similar 
colours and, hence, lead to the same level of agreement with 
the observed colour-colour distributions. There are many ob- 
servational evidences, from spectroscopic and photometric 
studies favouring the first g roup of solutio ns for the forma- 
tion of massive galaxies (e.g.[Renzini 200(|). Indeed, galaxies 
that formed a considerable amount of stars very recently, as 
in the second group, would show strong Hp absorption in 
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Figure Al. Best solution SFHs from the Montecarlo simulations described in Section 16.51 The first column shows the SFH in terms 
of fraction of mass formed at different ages. The second one shows the metallicity associated to each of the bursts, and the third one 
the slope of the IMF, assumed to be unimodal. In each panel the red, green and blue bars indicate the fractions corresponding to the 
contributions of the old, intermediate and young burst, respectively. The models are ordered according to the contribution of the young 
stellar populations. The vertical dashed line in each panel indicate, from left to right, an age of 5 Gyr, solar metallicity and fx=1.3 
(Salpeter IMF), respectively. We indicate in the last panel of each model the difference in u — g, g — r and r — i colours with respect to 
the median of the observed colours. 
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Figure Al. Continued 



